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ABSTRACT 


A  detailed  analysis  of  the  surface  wave  radiation  from 
two  underground  explosions  (BILBY  and  SHOAL)  and  an  earthquake 
(near  Fallon,  Nevada)  whose  epicenter  is  only  60  km  from  SHOAL 
indicates  that:  (1)  at  long  periods  the  surface  wave  radiation 
from  the  earthquake  can  be  explained  by  a  pure  quadrupole 
(double-couple)  source,  but  at  higher  frequencies  the  radiation 
pattern  shows  asymmetries  which  suggest  effects  due  to  rupture 
propagation  which  require  higher-order  multiple  terms  in  the 
source  equivalent  representation;  (2)  the  surface  waves  from 
the  explosions  can  be  explained  by  superposed  monopole  and 
quadrupole  sources,  with  no  indication  of  higher-order  multi- 
pole  terms,  at  least  in  the  frequency  range  comparable  to 
that  in  which  the  earthquake  signal  was  recorded;  (3)  a 
principal  conclusion  of  the  study  is  that  the  anomalous  radia¬ 
tion  from  explosions  is  probably  due  to  stress  relaxation 
around  the  shock-generated  shatter  zone  and  not  due 
to  earthquake  triggering. 

A  comparative  analysis  of  SHOAL  and  FALLON  shows  that:  (1) 
the  ratio  of  the  Love  wave  amplitude  generated  by  the  earthquake 
to  the  uove  wave  amplitude  from  the  explosion  is  proportional 
to  the  period,  which  implies  a  larger  source  dimension  for 
FALLON;  (2)  the  normalized  spectral  ratio  of  Love  wave  amplitude 


to  Rayleigh  wave  amplitude,  considered  as  a  function  of  period, 
is  near  unity  for  the  explosions  but  iarger  for  the  earth¬ 
quake  by  a  factor  of  two  or  three,  and  increasing  with  period. 
These  differences  might  be  useful  in  distinguishing  earth¬ 
quakes  from  explosions  (at  least  in  the  magnitude  range  of 
the  events  used  in  this  study,  4.4  and  above),  as  well  as 
for  estimating  source  parameters  (such  as  stress)  which  are 
of  fundamental  geophysical  interest 
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INTRODUCTION 


The  surface  waves  radiated  by  explosions  and  by  earthquakes 
arc  of  special  interest  since  they  provide  information  about 
the  intrinsic  character  of  the  source  and  the  propagation 
medium.  In  this  study  we  attempt  to  answer  the  following 
questions:  first,  what  are  the  spatial  radiation  patterns 
from  earthquakes  and  explosions,  and  how  do  they  differ? 

Second,  what  arc  the  differences  between  the  long-period  seismic 
radiation  from  these  two  types  of  sources,  as  determined 
from  the  surface  wave  spectra? 

These  questions  arc  related  to  the  problem  of  discrimi¬ 
nation  between  the  source  types,  and  arc  comparative  in 
nature.  In  addition  to  considering  these  questions,  we  also 
attempt  to  provide  a  theoretical  explanation  for  the  observed 
differences,  and  a  quantitative  description  of  the  individual 
source  properties. 

Since  this  is  primarily  a  comparative  study  intended  to 
show  differences  between  earthquakes  and  explosions,  it  is 
desirable  to  have  at  least  one  earthquake-explosion  pair  with 
a  common  origin  point  and  common  recording  stations,  so  that 
the  effects  of  propagation  can  be  minimized.  Ideally  the  sources 
should  also  be  matched  in  radiated  energy,  so  that  a  simple 
direct  comparison  can  be  made  without  the  necessity  for  correc¬ 
tions  which  arc  complicated  and  perhaps  not  completely  understood. 
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Of  course,  these  conditions  arc  impossible  in  practice;  the 
closest  we  could  come  to  meeting  them  was  to  use  records  from 
the  Long  Range  Seismic  Measurement  (l.RSM)  station  network, 
which  provides  good  azimuthal  and  distance  coverage  from 
the  Nevada  Test  Site,  with  well-calibrated  sensitive  instru¬ 
ments;  we  also  chose  an  earthquake  matching  one  of  the 
explosions  as  closely  as  possible  in  location  and  magnitude. 
We  found  that  an  appropriate  event  pair  for  this  purpose  is 
the  SHOAL  explosion  and  the  earthquake  of  20  .July  1962  at 
ballon,  Nevada.  The  parameters  of  all  the  events  used  in 
this  study  are  given  in  Table  1.  It  is  seen  that  the 
desirable  criteria  described  above  are  reasonably  well 
fulfilled* 

Since  we  arc  interested  in  the  nature  of  the  explosion 
source,  we  also  consider  the  IULBY  explosion,  which  was  well- 
recorded  and  of  reasonably  large  magnitude.  We  compare  this 
explosion  with  SHOAL  in  order  to  verify  that  the  observed 
properties  of  the  radiation  from  SHOAL  arc  typical  of  other 
explosion  events,  and  also  to  judge  the  variation  of  these 
properties  with  changing  yield  and  source  medium. 

figure  1  shows  the  location  of  the  sources  and  the 
stations.  It  is  evident  that  many  of  the  receivers  have  a 
common  propagation  path  for  the  SHOAL  and  l-ALLON  events. 

I igure  2  shows  the  response  of  the  long-period  instruments 


at  the  LRSM  stations.  The  frequency  range  in  which  the  signal- 
to-noisc  ratios  required  for  this  research  were  deemed  adequ.it 
is  indicated  by  the  horizontal  bar;  all  spectral  data  were 
corrected  for  instrument  response  within  this  range,  and 
data  outside  this  range  were  rejected. 

Unfortunately,  the  azimuthal  coverage  of  stations  toward 
the  west  of  the  sources  is  rather  sparse,  and  therefore  some 
of  our  fits  to  the  azimuth-dependent  data  must  be  regarded  as 
tentative.  Nevertheless,  a  large  amount  of  data  is  available, 
and  we  can  support  our  specific  (albeit  tentative)  conclusions 
by  a  reasonably  good  general  agreement  with  our  theoretical 
predictions . 

Part  of  the  theoretical  framework  of  this  study  is  based 
on  the  work  of  Archambcau  (1968;  1971).  Briefly,  this  theory 
provides  a  description  of  tectonic  effects  in  terms  of  an 
initial-value  formulation,  taking  into  account  the  initial 
stress  field  and  its  relaxation  during  the  rupture  process. 
Sources  of  this  type  are  therefore  termed  relaxation  sources. 

A  discussion  of  the  theory  in  the  context  of  the  present 
study  is  given  by  Archambeau  (1971)  in  greater  detail 
than  we  need  present  here.  The  present  application  of  the 
theory  is  two-fold:  first,  as  a  means  of  guiding  the  develop¬ 
ment  of  this  study  and  in  interpreting  the  observations  from 
the  earthquake  source;  second,  as  a  framework  for  the 


interpretation  of 
sources,  which  we 
later  section  the 
will  be  discussed 


the  anomalous  radiation  from  explosive 
consider  to  be  of  tectonic  origin.  In  a 
relevance  and  application  of  this  theory 
in  greater  detail. 
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SIGNAL  ANALYSIS 


tVe  have  studied  the  surface  wave  radiation  from  the 
three  seismic  events,  using  a  digital  computer  program 
which  measures  both  the  amplitude  spectrum  and  the  dispersion 
of  the  surface  waves.  This  program  is  based  on  programs  origi¬ 
nally  written  by  Alexander  (11)63),  and  is  described  in  detail  l>\ 
Arehanbcau  and  I  linn  (1965).  The  program  automatically  performs  ,i 
sequence  of  operations,  ultimately  giving  the  signal  group 
velocity  and  amplitude  spectrum.  The  amplitude  spectrum  is 
compared  to  a  sample  of  noise  preceding  or  following  the 
signal,  and  the  program  flags  and  deletes  from  further 
analysis  the  period  ranges  in  which  the  signal-to-noise 
ratio  is  too  small  to  allow  meaningful  interpretation  of 
the  signal  spectrum.  The  group  velocity  is  calculated  using 
a  comb  of  narrow-hand  recursive  digital  filters  (Archambeau 
and  FI inn ,  1965).  At  each  of  many  closely  spaced  frequencies 
the  signal  data  window  is  passed  through  a  narrow-band  phase- 
less  digital  filter  (Q  =  500  in  this  study;  for  background  see 
also  Rader  and  Gold,  1969).  The  envelope  of  the  filter  output 
at  each  frequency  is  calculated,  and  the  time  of  occurrence 
of  the  largest  maximum  in  the  envelope  is  picked.  This 
maximum  occurs  at  the  group  arrival  time  (i.e.,  the  energy 
arrival  time)  of  the  dominant  mode  present  at  that  frequency. 

The  cpiccntral  distance  divided  by  this  time  is  thus  the 
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ipparent  group  velocity  of  the  dominant  mode,  at  the  center 
frequency  of  the  filter.  Group  velocities  of  secondary 
envelope  maxima  are  also  detected  and  displayed. 

The  variation  with  frequency  of  the  amplitude  at  the 
time  of  occurrence  of  the  dominant  mode  arrival  gives  an 
estimate  of  t lie  amplitude  spectrum  of  the  dominant  mode, 
and  this  therefore  provides  a  method  for  determining 
surface  wave  mode  spectra  as  a  function  of  frequency.  We 
call  such  spectra  group  spectra.  These  can  be  determined 
for  cacli  mode  present  at  a  given  frequency,  Figure  3  shows 
examples  of  raw  data  and  the  group  velocity  curves  determined 
hv  the  program;  figure  4  shows  a  comparison  between  the 
ordinary  Fourier  amplitude  spectra  and  the  group  spectra. 

We  have  found  in  general  that  the  two  estimates  agree  to 
within  a  few  percent  in  the  frequency  range  in  which  the 
s i gnal-to-noisc  ratio  is  greater  than  1,5,  even  when  more 
than  one  mode  is  present.  We  believe  the  group  spectrum  to 
he  intrinsically  more  reliable  than  the  Fourier  spectrum, 
since  we  can  associate  a  group  arrival  time  with  each 
spectral  estimate,  and  thus  eliminate  those  points  in 
the  spectrum  which  are  badly  contaminated  by  noise.  We 
have  therefore  used  the  group  spectral  estimates  consistently 
in  the  present  study. 
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To  calculate  the  theoretical  Love  and  Rayleigh  spectra, 
it  is  necessary  to  make  some  assumption  about  the  veloc  ity 
structure  in  the  vicinity  of  the  source.  For  the  structure  near 
BILBY  we  used  model  35CM2  of  Alexander  (1963).  For  SHOAL  and 
Fallon  we  used  the  structure  previously  used  by  Toksoz  ct  al. 
(1965),  which  is  a  modification  of  baton's  (1963)  structure 
between  San  Francisco  and  liurcka,  Nevada. 
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CONCEPTS  AND  THEORY 


An  important  concept  underlying  the  theoretical  frame¬ 
work  of  this  research  is  that  the  radiated  elastodynamic 
field,  outside  the  region  containing  the  source,  can  be 
represented  in  the  frequency  domain  by  a  superposition  of 
multipole  contributions.  This  fact,  familiar  in  the  analogous 
cases  of  potential  theory  and  electrodynamics,  was  formulated 
by  Archambeau  (1968)  for  an  elastic  radiation  field.  The 
implication  is  that  we  can  represent  the  source  of  the 
radiated  field  as  an  equivalent  point  source  located  anywhere 
inside  the  actual  region  of  energy  release,  and  we  can 
adjust  the  multipole  coefficients  so  that  the  observed  field 
outside  the  source  region  is  identical  to  that  caused  by 
the  true  source.  Since  the  multipoles  are  individually  equiva¬ 
lent  to  particular  superpositions  of  point  forces,  we  can  (if 
we  wish)  speak  of  force  couples,  double  couples,  etc.,  instead 
of  multipoles,  quadrupoles,  and  so  forth.  In  the  present  study 
we  discuss  the  equivalent  point  source,  for  either  the 
explosions  or  the  earthquakes  in  terms  both  of  the  appropriate 
multipoles  and  their  point  force  equivalents,  depending  on 
which  designation  is  more  familiar  and  useful  in  the  particular 
context . 
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For  the  purpose  of  this  report  we  will  merely  summarize 
the  properties  of  the  theory  for  the  specific  cases  of  interest, 
and  then  use  these  properties  in  the  particular  applications 
we  wish  to  make.  More  detail  on  the  theory  is  given  in  Appendix  I 

I.  For  the  shock-induced  spherical  shatter  zone,  the  tec¬ 
tonic  field  is  a  pure  quadrupole  (double  couple)  field,  and  the 
frequency  dependence  of  this  quadrupole  field  is  similar  to  that 
of  the  monopole  component  caused  by  the  explosion  itself,  since 
the  source  dimensions  are  similar.  The  spectra  for  both  source 
components  are  both  quite  sharply  peaked  at  nearly  the  same 
frequency,  and  the  spectra  fall  off  with  similar  shapes  on  both 
sides  of  the  peak  (see  Archambeau,  1971,  for  further  details). 
These  general  results  hold  both  when  the  rupture  velocity  ir 
less  than  or  equal  to  the  shear  velocity. 

II.  For  the  propagating  spherical  rupture,  which  we  adopt 

as  our  model  of  the  earthquake  source  at  frequencies  sufficiently 
low  that  the  wave  length  of  the  radiated  wave  field  is  consider¬ 
ably  larger  than  the  rupture  dimensions,  the  quadrupole  or  double 
couple  terms  heavily  dominate  the  solution,  and  the  field  is 
essentially  pure  quadrupole.  At  higher  frequencies  the  higher- 
order  multipoles  can  become  comparable  in  amplitude.  This  is  a 
result  of  the  rupture  propagation  effect,  and  these  multipole 
contributions  are  such  that  the  spectrum  has  a  sequence  of 
maxima  and  minima  at  high  frequencies,  and  the  spatial  radiation 
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pattern  is  more  complex  than  that  for  a  pure  quadrupole  pattern. 
In  particular,  the  pattern  has  the  general  form  of  a  quadrupole 
pattern  but  with  larger  amplitudes  either  in  the  direction  of 
rupture  propagation  or  in  the  opposite  direction  (see  Archambeau 
and  Minster,  1971,  for  details).  The  spectrum  of  surface  waves 
from  this  type  of  source  has  an  absolute  maximum,  the  location 
and  amplitude  of  which  is  controlled  primarily  by  the  length 
of  the  zone  swept  out  by  the  propagating  rupture.  The  average 
falloff  of  the  spectrum  on  both  sides  of  this  peak  is  not  as 
steep  as  fov  either  the  individual  monopole  or  quadrupole 
componets  of  the  spectrum  of  surface  waves  from  an  explosion. 

In  view  of  these  results  we  will  use  the  following  approach 
in  the  interpretation  of  the  explosion  and  earthquake  radiation 
fields  : 

(1)  An  explosion  in  a  prestressed  medium  will  be  repre¬ 
sented  by  superposed  monopole  and  quadrupole  fields.  Since  the 
quadrupole  field  is  equivalent  to  a  double-couple  point  force 
system,  we  will  adopt  the  force  equivalent  convention  in 
discussing  the  tectonic  component  of  this  source.  We  see  from 
our  previous  discussion  that  if  the  higher-order  multipoles 
.re  of  tectonic  origin,  the  quadrupole  component  will  be 
essentially  the  only  term  worth  considering  for  Love  waves  at 
long  periods.  Since  we  will  analyze  the  surface  wave  radiation 
in  the  "long"  period  range,  then  whatever  type  of  tectonic 
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effect  is  involved  --  whether  triggering  of  a  linear  rupture 
or  simple  relaxation  around  the  shock- induced  shatter  zone  --  the 
quadrupole  field  will  still  be  dominant. 

If  induced  faulting  on  a  large  scale  does  occur,  then  we 
would  expect  to  see  rupture  propagation  effects  at  the  short- 
period  end  of  our  spectra.  Further,  if  such  triggering  occurs, 
the  pattern  shape  for  Rayleigh  waves  at  long  periods  should 
appear  to  be  nearly  an  ideal  quadrupole  pattern  rather  than  a 
superposed  quadrupole-monopole  pattern,  since  the  quadrupole 
excitation  for  large-scale  faulting  will  be  considerably  greater 
at  very  long  periods  than  will  the  monopole  excitation  due  to 
the  explosion,  simply  because  of  the  difference  in  source 
dimensions  for  the  two  source  types.  If  the  anomalous  component 
of  the  field  is  associated  only  with  relaxation  around  a  roughly 
spherical  shatter  zone,  the  Rayleigh  wave  radiation  pattern 
should  have  nearly  the  same  shape  over  the  entire  frequency 
range  of  observation  because  of  the  similarity  in  multipole 
spectral  excitation  for  the  two  source  types;  and  the  entirely 
anomalous  Love  wave  pattern  should  not  show  rupture  propagation 
effects:  it  should  be  pure  quadrupole.  We  will  examine  the 
spectral  observations  from  two  explosions  with  a  view  to 
ascertaining  which  of  these  mechanisms  is  operative  in  a  given 
case.  Of  course,  we  will  also  entertain  the  possiblity  that  the 
anomalous  field  is  not  tectonic  at  all,  but  is  conversion 
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due  to  cracking  or  non-parallel  layering.  In  this  study  our 
main  argument  for  acceptance  of  tectonic  release  rather  than 
conversion  is  the  quality  of  the  fit  to  the  observations  that 
this  hypothesis  furnishes.  At  the  very  least,  the  observational 
data  can  be  taken  at  face  value,  and  we  can  argue  that  we  have 
quantitatively  determined  the  first  two  multipole  contributions, 
whatever  their  origin  might  be. 

(2)  The  quadrupole  component  of  the  surface  wave  radiation 
from  the  earthquake  will  be  studied  quantitatively.  Hence  we 
represent  the  earthquake  by  a  double-couple  point  source,  and 
us  we  have  already  seen,  this  will  be  the  appropriate  representa¬ 
tion  at  long  periods.  At  the  short-period  end  of  the  surface 
wave  spectrum  we  may  anticipate  seeing  higher-order  multipolc 
effects  due  to  rupture  propagation.  Since  our  range  of  observa¬ 
tions  is  at  rather  long  periods,  the  double-couple  source  will 
turn  out  to  provide  an  appropriate  spatial  radiation  pattern 
which  can  be  used  to  fit  the  observations.  Depending  on  the 
rupture  length,  however,  we  might  expect  to  see  departures  in 
the  radiation  pattern  from  the  pure  quadrupole  type,  as  well 
as  the  series  of  maxima  and  minima  toward  the  short-period  end 
of  the  spectrum,  due  to  rupture  propagation.  We  shall  there¬ 
fore  check  our  observations  with  this  in  mind,  although  in  the 
present  study  we  make  no  attempt  to  fit  the  details  of  the 
source  spectrum,  nor  try  to  study  the  effect  of  higher-order 
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multipole  terms  in  the  radiation  patterns.  Rather,  we  simply 
fit  the  spatial  radiation  pattern  up  to  the  quadrupolc  term, 
and  draw  conclusions  from  the  nature  of  the  fit.  We  then 
qualitatively  compare  the  details  of  tltc  observed  spectrum 
and  the  observed  spatial  radiation  patterns  to  the  theoretical 
predictions. 

Since  we  will  deal  in  detail  with  monopolc  and  quadrupo le 
fields  in  this  study,  it  is  convenient  to  introduce  parameters 
related  to  their  force  equivalents.  For  the  double  couple  (the 
quadrupolc  equivalent)  we  define  a  plane  along  which  one  of 
the  couples  is  oriented,  the  other  couple  being  oriented  in  a 
direction  perpendicular  to  it.  Figure  5  shows  such  a  plane,  wi 
its  orientation  in  space  specified  by  the  unit  vectors  shown 
For  an  earthquake,  this  plane  corresponds  to  the  fault  plane, 
which  is  simp  .y  a  limiting  form  of  the  surface  enclosing  the 
region  of  failure.  If  this  region  has  a  small  width  relative 
to  its  other  dimensions,  there  is  little  geometrical  distinc¬ 
tion  to  be  made  between  the  envelope  surface  and  the  fault 
plane  shown  in  Figure  5. 

As  is  customary,  the  orientation  of  the  fault  plane  is 
described  by  a  strike  angle,  measured  from  north,  and  a  dip 
angle  measured  from  the  normal  to  the  strike  direction,  as 
shown.  A  point  of  observation  (0)  will  he  specified  by  an 
azimuth  angle  0,  measured  from  the  strike  direction,  and  an 


epicentral  distance,  r,  measured  from  the  point  of  initiation  of 
t lie  rupture;  this  latter  point  is  shown  in  Figure  5  as  the 
origin  of  coordinates.  If  the  fault  plane  or  actual  rupture 
surface  nay  terminate  at  some  depth  h  below  the  earth's  surface, 
this  distance  is  measured  in  the  direction  of  axis  x^  in  Figure  5, 
normal  to  the  surface  of  the  earth. 

he  define  t/*^  in  Figure  5  as  the  traction  vector  at  a 
distance  -c/2  just  behind  the  fault  plane,  when  the  fault 
plane  is  viewed  from  the  direction  defined  by  n.  When  the  fault 

A 

plane  is  approached  from  the  -n  direction  there  is  an  equal 
and  opposite  traction  -t^0^  at  a  point  a  distance  ♦c/2  just  in 
front  of  the  plane.  This  configuration  of  tractions  defines  one 
of  the  force  couples  for  the  equivalent  source;  its  orientation 
in  the  plane  will  be  described  by  the  so-called  slip  angle, 
which  is  simply  the  angle  \  defined  with  reference  to  t ^ ,  as 
shown  in  Figure  5.  For  the  full  quadrupole  equivalent  we  require 
an  equal  and  opposite  couple  whose  orientation  is  perpendicular 
to  the  plane  shown  in  Figure  5,  such  that  the  net  torque  at  the 
point  I*  is  zero.  This  particular  configuration  is  required  by 
the  equilibrium  equations,  since  for  a  relaxation  process  no 
net  torques  can  be  maintained.  The  origin  of  the  second  couple 
is  the  Poisson  effect  which  occurs  during  relaxation. 

In  the  ease  of  an  earthquake  the  definition  of  the  fault 
plane  and  the  system  of  couples  just  described  is  intuitively 
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appealing  and  reasonably  straightforward  (although  from  a  purely 
intuitive  approach,  only  a  single  couple  night  appear  to  be 
required  at  first  sight;  this  fact  has  contributed  to  the  well- 
documented  disagreement  among  seismologists  as  to  whether  a 
single  or  double  couple  source  is  required).  For  stress  relaxa¬ 
tion  associated  with  a  "spherical"  shatter  tone,  however,  the 
same  equivalent  double-couple  point  source  applies,  but  of 
course  the  physical  identification  of  a  "fault  plane”  is  no 
longer  obvious.  Even  so,  we  will  define  a  plane  through  the 
source  and  its  associated  orientation  parameters  for  this 
problem  as  well.  In  this  case,  which  may  apply  to  explosion¬ 
generated  earthquakes  in  the  real  world,  the  plane  and  the 
couple  orientations  are  related  to  the  orientation  of  the 
prestress  field. 

This  approach  allows  us  to  use  Ben-Menahem  and  Harkrider's 
(1964)  point  source  results  directly,  to  describe  the  quadru- 
pole  field  from  an  earthquake  and  from  an  explosion  in  a 
stressed  medium.  Since  we  are  primarily  interested  in  theo¬ 
retical  fits  to  the  azimuthal  dependence  of  the  radiation  fields, 
we  need  only  calculate  the  variation  of  a  single  function,  x(°) 
which  gives  the  radiation  patterns  for  either  Love  or  Rayleigh 
waves  from  a  double-couple  point  source  at  a  depth  h  in  a 
plane-layered  earth  model.  In  particular,  Ben-Menahem  and 
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Ilarkrider  ( 1 9  <»  4  ;  equation  25)  give: 


X ( 0 )  a  J0  *  i  (d A s in  0  ♦  d2cos  0)  ♦  dj  sin  20  ♦  d4  cos  26 

where  for  Rayleigh  waves  the  coefficients  are: 

dQ  *  1/2  B(h)  sinX  sin  26 

d^  =  -  C ( h )  sinX  cos  26 

d^  *  -  C ( h)  cosX  cos  6 

d,  =  A (h)  cosX  sin  6 

d4  ■  -1/2  A (h)  sinX  sin  26 

and  for  Love  waves  the  coefficients  are: 

d„  -  0 
o 

dj  ■  G(h)  cos  X  cos  6 
d^  =  -G(h)  sin  X  cos  26 
dj  ■  1/2  sin  X  sin  26 
d4  »  cos  X  sin  6 

The  functions  A(h),  B(h),  C(h),  and  G(h)  depend  on  the 
equivalent  point  source  depth  and  on  the  details  of  the 
earth  model  chosen.  They  are  given  explicitly  by  Ben-Menahem 
and  Harkridcr  (1964;  equations  26-36),  but  since  they  are  quite 
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lengthy,  they  will  not  he  reproduced  here.  We  will  vary  , 

6,  and  li  to  lit  our  ohscrv.it  i  oils  ,  while  in  a  i  nta  i  11  i  n>  i  I  i.u  1 
average  crust-upper  mantle  model  appropriate  to  the  wstnn 
United  States  i  a  the  remainder  of  this  work.  Without  detail*.  1 
path  corrections  (or  control),  therefore,  we  expect  to  see 
scatter  m  the  observations  about  any  fit  we  might  achieve, 
this  will  to  some  extent  limit  our  ability  to  draw  hard  conclu¬ 
sions.  However,  with  a  good  average  structural  model  we  can 
expect  to  resolve  and  identify  all  hut  the  more  subtle  radia¬ 
tion  pattern  variations. 

For  the  monopoie  field  from  the  explosion  we  merely  not< 
that  it  has  the  form  (in  our  present  notation): 

x'2)d.  “>  ■  Ou>  "o2Xt' 


with 


Au  (u) 
oo v  ' 


(I 

(l 

0 

fi(u) 


The  function  (J(m)  will  depend  on  the  radius  of  the  elastic 
plastic  boundary  and  the  effective  pressure  at  this  radius. 
Archamheau  (1D7I)  describes  such  functions  and  the  various 
ways  of  determining  them.  It  is  clear,  however,  that  we  can 
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again  use  a  point  source  model  and  employ  Ben-Mcnahcm  and 
ilarkrider's  results  for  a  point  pressure  source  without 
being  directly  concerned  with  all  the  details  of  the  function 
f.(  ) . 

'.n  [.articular,  our  procedure  in  dealing  with  the  explosion 
is  to  superpose  a  point  pressure  and  point  double-couple  sources 
and  to  scale  the  double-couple  component  by  a  factor  I  (see 
Toksoz  et  al.,  lt)f>5).  lor  observations  of  the  radiation  pattern 
at  a  given  frequency,  then,  we  can  adjust  this  factor  1,  along 
with  the  other  parameters  of  the  double  couple,  in  order  to 
obtain  a  theoretical  fit  to  the  observed  radiation  pattern.  If 
ti ( («: )  is  such  that  the  actual  source  monopole  spectrum  is 
similar  to  the  quadrupole  spectrum  over  the  range  of  observa¬ 
tion,  then  the  fit  achieved  at  one  frequency  should  hold  over 
the  entire  frequency  range.  If  the  monopole  excitation  is  quite 
different  from  that  associated  with  the  quadrupole  term,  then 
the  1  factor  will  change  with  frequency.  From  our  earlier 
discussion  of  the  effects  cif  shock- induced  relaxation  we  note 
that  if  relaxation  around  the  shatter  zone  is  all  that  is 
involved  tectonically,  then  the  F  factor,  once  determined, 
should  he  found  to  remain  fairly  constant  with  frequency. 

If  triggering  of  large-scale  asymmetric  faulting  is  involved, 
the  F  factor  should  increase  strongly  with  increasing  period, 
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reflecting  greater  quadrupole  excitation  at  long  periods. 

Thus  we  will  attempt  to  fit  radiation  patterns  over  a  range 
of  frequency  in  order  to  determine  whether  such  variations 
in  F  exist* 

To  conclude  this  part  of  our  discussion,  it  is  worth¬ 
while  to  note  some  of  the  less  obvious  difficulties  inherent 
in  fitting  theoretical  predictions  to  the  observations. 

These  difficulties  arise  primarily  from  a  lack  of  uniqueness 
in  fitting  the  observations,  and  are  best  illustrated  by  a 
few  examples. 

In  particular,  since  the  azimuthal  variation  of  the 
ratio  of  Love  to  Rayleigh  wave  amplitudes  (L/R'  n  a  given 
frequency  is  an  important  variable,  it  is  appropriate  at 
this  point  to  discuss  the  effect  on  the  L/R  radiation  pattern 
of  variations  in  the  free  parameters  describing  the  source. 

For  an  explosion  with  an  anomalous  quadrupole  or  double¬ 
couple  component,  we  have  as  variables  the  depth  of  the 
double  couple;  the  strike,  dip,  and  slip  angles  describing 
the  orientation  of  the  double  couple;  and  the  F  factor 
relating  the  monopole  excitation  to  the  quadrupole  excitation. 

We  note  parenthetically  that  the  depth  of  the  equivalent 
double  couple  required  to  explain  the  observations  could, 

*We  note  also  that  intuition  suggests  a  roughly  constant  F  value 
if  conversion  by  cracking  causes  the  Love  waves. 
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if  it  is  significantly  different  from  the  known  depth  of  the 
explosions,  be  diagnostic  of  the  mechanism  of  the  quadrupole 
field:  a  much  greater  depth,  for  example,  would  suggest 

triggering  of  an  earthquake,  while  a  depth  near  that  of  the 
monopole  equivalent  might  be  caused  either  by  triggering  or 
by  relaxation  due  to  the  shatter  zone  itself  or  by  conversions. 

For  an  earthquake  we  have  the  same  parameters  available 
except  for  the  F  factor,  since  for  earthquakes  we  assume 
that  there  is  no  monopole  contribution. 

Many  variations  of  combinations  of  these  parameters  can 
give  rise  to  closely  similar  L/R  patterns.  Figure  6  shows, 
for  example,  the  variation  in  L/R  pattern  for  an  explosion 
in  a  stressed  medium  for  different  double -couple  slip  angles. 

A  similar  variation  in  the  pattern,  involving  slight 
azimuthal  rotations  and  changes  of  magnitude,  can  be  caused 
by  changes  in  the  depth  of  the  double-couple  component. 
Similarly,  radically  different  values  of  slip  angle  and 
double-couple  depth  equivalent  can  result  in  the  same 
observed  radiation  pattern.  A  similar  ambiguity  exists  in 
the  combination  of  dip  angle  and  the  F  factor  for  an  explosion. 
Figure  7  shows  the  variation  In  the  L/R  pattern  for  varying 
dip  angle,  holding  the  other  parameters  fixed.  Here  the 
pattern  lobes  increase  in  magnitude  with  increasing  dip,  and 
essentially  the  same  kind  of  variation  occurs  as  F  increases. 
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Furthermore,  various  combinations  of  strike,  slip,  and  dip 
angles,  and  the  F  factor  can  be  found  that  give  very  similar 
radiation  patterns,  so  that  we  are  confronted  with  the 
fundamental  lack  of  uniqueness  which  is  characteristic  of  most 
geophysical  problems.  In  our  treatment  of  the  actual 
observations  and  the  theoretical  fit  to  our  data,  we  there¬ 
fore  work  in  a  period  range  in  which  the  scatter  of  data  is 
small,  so  that  we  can  with  reasonable  justification  make 
relatively  fine  distinctions  in  the  fitted  models.  We  also 
constrain  our  solutions  to  be  in  at  least  rough  agreement 
with  the  known  tectonic  setting  of  the  source  area:  thus 
we  make  use  of  the  fact  that  something  is  known  of  the 
strikes  and  dips  of  the  faults  in  the  source  region,  and  wc 
impose  a  set  of  constraints  on  our  allowable  source  models. 
Further,  the  ambiguities  for  the  earthquake  are  not  as  serious 
as  for  the  explosions,  since  wc  deal  with  only  one  type  of 
equivalent  source,  and  wc  use  the  fit  for  the  double  couple 
of  the  earthquake  as  a  guide  for  possible  solutions  for  the 
SHOAL  explosion.  This  presupposes  tectonic  release  as  the 
origin  of  the  double -couple  component,  of  course,  and  hence 
only  serves  to  constrain  permissible  solutions  of  this  type. 
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SURF  ACL  WAV  I:  RADIATION  PATTERNS  AND  DEDUCTION  OF  THE 
EQUIVALENT  MULTI  POLE  SOURCE 

In  this  section  we  consider  the  surface  wave  radiation 
patterns  associated  with  the  two  explosions  and  the  single 
earthquake  studied,  and  discuss  the  point-source  models 
implied  by  the  data. 

Figure  8  shows  Love  and  Rayleigh  wave  patterns  for  two 
BILBY  explosion  frequencies  spanning  the  bandwidth  in  which 
the  signal  power  was  greatest.  We  note  that  the  s tati ons  closest 
to  the  sources  give  an  adequate  azimuthal  coverage  for  our 
purposes,  although  more  stations  to  the  west  would  have  been 
desirable.  It  is  reasonably  clear  that  variations  in  crust 
and  upper  mantle  structure  strongly  perturb  the  radiation 
patterns;  in  particular,  the  stations  in  the  eastern  United 
States  show  quite  large  amplitudes  which  consequently  cause 
distortions  in  the  contours  of  the  pattern.  It  appears  that 
eastward-propagating  surface  wave  energy,  once  across  the 
Rocky  Mountains,  propagates  from  there  on  with  little 
attenuation,  and  that  the  observed  surface  wave  amplitudes 
can  actually  be  greater  in  the  East  than  at  stations  in  the 
western  tectonic  provinces  which  are  much  closer  to  the  source. 
This  can  be  explained  in  part  by  the  fact  that  the  crust  in 
tlie  eastern  United  States  has  intrinsically  higher  wave 


-22- 


velocities  in  the  intermediate  and  deep  layers  than  does  the 
crust  in  the  Basin  and  Range  Province  (see,  for  example, 
Archambeau  ct  al . ,  1969) .  The  top  of  the  upper  mantle  also 
has  higher  velocity  in  the  East  than  in  the  West.  This  results 
in  a  distribution  of  the  surface  wave  energy  nearer  the  surface 
for  surface  waves  in  the  period  range  observed  (5  to  10  seconds) , 
and  hence  hither  amplitudes  at  the  surface,  even  though  the 
total  wave  energy  is  less.  Further,  the  lateral  variations  within 
the  eastern  provinces  appear  to  be  milder  than  those  in  the 
tectonic  provinces  in  th  West  (Archambeau  et  al.,  1968,  1969) 
and  hence  there  is  less  scattering  and  more  uniform  propagation 
of  the  waves  generally. 

The  patterns  in  Figure  8  for  both  Love  and  Rayleigh  waves 
can  be  approximated  by  superposed  quadrupole-monopole  sources 
with  both  equivalent  point  sources  at  the  depth  of  the  detona¬ 
tion  point  for  the  explosion.  These  theoretical  patterns, 
approximately  the  same  for  T  =  15  or  20,  are  shown  in  the  insets 
in  Figure  8.  The  source  parameters  are  as  given  in  the  figure 
and  listed  in  Table  II;  the  ellipticity  parameter  e  is  deter¬ 
mined  by  the  structure  used  to  predict  the  radiation  patterns. 

We  were  particularly  interested  in  the  effect  of  varying  the 
depth  of  the  double  couple  equivalent,  but  found  that  the  best 
fit  to  the  observations  is  given  by  a  double  couple  at  the 
depth  of  the  explosion.  This  supports  a  local  effect  such  as 
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tectonic  release  or  conversion  in  ths  very  near  vicinity  of 
the  explosion,  but  does  not  clearly  distinguish  between  induced 
faulting  as  opposed  to  relaxation  due  only  to  the  shock-generated 
"spherical"  fracture  zone. 

It  is  also  significant  that  the  same  combination  of  mono¬ 
pole  and  quadrupole  sources  fits  both  the  Love  and  Rayleigh  waves. 
Further,  we  found  very  large  changes  in  observed  amplitudes  of 
the  radiation  patterns  between  two  frequencies  where  the  shape 
of  the  pattern  remains  essentially  invariant.  In  addition,  the 
pattern  did  not  vary  significantly  from  that  shown  over  the 
entire  range  from  approximately  10  seconds  to  35  seconds,  although 
at  the  extremes  of  this  range  the  observed  power  was  low  and  the 
pattern  resolution  became  poor.  Nevertheless,  we  conclude  that 
there  is  no  evidence  of  pattern  shape  change  over  this  range. 

This  supports  the  idea  that  the  mechanism  associated  with  the 
double  couple  source  is  not  only  located  in  the  near  vicinity  of 
the  explosion,  but  that  the  spectrum  produced  by  this  mechanism 
is  similar  in  shape  (as  a  function  of  frequency  over  the  range 
observed)  to  the  explosion  spectrum.  Indeed,  the  scaling  factor 
F  =  0.5  for  the  double  couple  is  independent  of  frequency  over 
the  observed  range,  as  is  predicted  for  the  radiation  due  to 
relaxation  around  the  shock-created  shatter  zone.  The  magnitude 
of  this  excitation  factor  is  also  of  the  order  expected  for 
relaxation. 
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The  amplitude  and  energy  spectrum  produced  by  this  double 
couple  source  component  from  Bilby  has  been  studied  by  Archambeau 
and  Sammis  (1970)  under  the  hypothesis  that  this  contribution 
to  the  field  is  actually  due  to  tectonic  prestress  relaxation 
around  the  shock- induced  shatter  zone.  It  was  found  that  a 
fit  could  be  achieved  for  a  homogeneous  shear  prestress  of 
70  bars  and  a  shatter  zone  of  roughly  400  meters  in  radius, 
which  are  reasonable  values.  Thus  we  tentatively  conclude  that 
the  observations  from  BILBY  support  the  hypothesis  for  this 
particular  form  of  tectonic  release.  We  do  not  see  a  change  in 
\  attern  shape  as  predicted  for  large-scale  faulting,  and  while 
we  cannot  definitely  rule  out  a  triggering  effect,  there  is 
little  evidence  to  support  it.  Aki  et  al.  (1969)  have,  however, 
drawn  the  opposite  conclusion  based  on  a  different  set  of 
observations . 

As  to  the  fundamental  question  of  whether  tectonic  effects 
are  responsible  for  the  double  couple  component,  we  note  at  this 
point  that  the  strike  and  dip  of  the  "couple  plane"  ($  =  342°, 

6  ■  90°)  are  consistent  with  the  orientation  of  tectonic  features 
in  the  vicinity  of  the  explosion  site.  This  simply  indicates 
that  the  average  stress  field  in  the  region  is  oriented  so  as  to 
produce  a  double  couple  field  in  agreement  with  the  observations 
for  either  fault  triggering  or  shatter-zone  relaxation  effects  -- 
therefore,  the  tectonic  hypothesis  is  at  the  very  least  consistent 
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with  our  observations.  Furthermore,  while  our  amplitude 
pattern  fit  indicates  a  double  couple  source  in  the  near 
vicinity  of  the  explosion  (within  a  few  kilometers),  observa¬ 
tions  of  the  Love  wave  group  velocity  requires  that  the  double 
couple  be  very  close  to  the  explosion  site,  in  order  that  the 
observed  group  velocities  be  compatible  with  our  rough  knowledge 
of  the  crustal  structure  along  the  propagation  paths.  Hence, 
while  one  might  argue  from  the  observations  of  the  Rayleigh 
waves  alone  that  the  patterns  might  be  accounted  for  by  ampli¬ 
tude  fluctuations  caused  by  local  structural  variations  at  the 
individual  stations,  it  is  impossible  under  such  a  hypothesis 
to  explain  the  Love  and  Rayleigh  wave  group  velocity  observa¬ 
tions,  which  arc  consistent  with  our  prior  knowledge  of  source- 
to-station  velocity  structure  within  the  Basin  and  Range 
Province.  In  particular,  only  a  double  couple  source  at  the 
explosion  location  can  explain  all  our  observations. 

The  theoretical  patterns  shown  in  Figure  8  provide  a  rather 
qualitative  way  of  fitting  the  observations;  that  is,  we  have 
tried  to  fit  by  comparing  the  general  shapes  with  the  contours 
obtained  from  the  observations.  A  more  quantitative  approach 
is  achieved  by  the  method  illustrated  in  Figure  9.  In  this  case 
we  can  simultaneously  consider  Love  and  Rayleigh  wave  amplitude 
patterns  (the  vertical  component  amplitude  in  the  case  of 
Rayleigh  waves)  by  fitting  their  ratio,  which  we  will  designate 


-26 


as  L/R.  This  can  be  done  as  a  function  of  frequency,  so  that  we 
consider  a  reasonable  frequency  sampling  of  such  L/R  patterns 
and  fit  them  with  a  fixed  combination  of  monopole  and  quadru- 
pole  parameters.  The  pawmeters  are  the  same  as  for  the  patterns 
shown  in  Figure  8  and  listed  in  Table  II,  of  course,  rnd  we 
only  note  that  the  L/R  fit  can  be  seen  to  be  quite  good.  All 
the  source  parameters  listed  remained  nearly  constant  over  the 
range  10  to  35  seconds.  We  will  return  to  the  spectral  character 
of  the  L/R  ratio  in  the  next  section,  but  at  this  point  it  is 
reasonable  to  strongly  suspect  that  the  mechanism  producing 
the  double  couple  component  for  BILBY  is  tectonic  stress  relaxa¬ 
tion  near  the  roughly  spherical  shatter  zone,  because  of  this 
constancy  of  source  parameters  versus  frequency. 

The  SHOAL  explosion  was  considerably  smaller  than  the  BILBY 
event  (Table  I),  but  a  double-couple  component  similar  to  that 
obtained  for  BILBY  was  observed.  In  particular,  Figure  10  shows 
radiation  patterns  similar  in  shape  to  those  observed  for  BILBY 
and  the  theoretical  fit  achieved  gave  an  F  factor  and  double 
couple  source  depth  which  are  like  those  for  BILBY.  The  same 
inteipretation  follows;  that  is,  the  stability,  or  shape 
invariance,  of  the  patterns  with  frequency  (over  the  range  10  to 
30  seconds  in  this  case) ,  the  coincidence  of  the  quadrupole  and 
monopole  source  depths,  the  strike  and  dip  of  the  equivalent 
"fault"  plane  coinciding  with  tectonic  features  of  the  area  and 
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the  moderate  constant  value  for  F,  all  imply  tectonic  relase  due 
to  stress  relaxation  near  the  explosive  shatter  zone.  Further, 
the  source  parameters  defining  the  double  couple  component  are 
nearly  identical  with  those  obtained  for  Bilby,  which  implies 
a  similar  tectonic  stress  field  in  the  two  areas? 

Figure  11  shows  the  L/R  pattern  at  a  period  of  15  seconds. 
The  fit  to  the  data  is  quite  good  considering  the  large 
variance  to  be  expected  in  the  ratio  of  two  spectra,  due,  for 
example,  to  structural  effects  and  noise.  The  scatter  from  the 
theoretical  fit  would  also  manifest  itself  in  the  form  of  dis¬ 
placements  of  the  observed  points  on  this  plot  in  azimuth,  due 
to  refraction  of  the  surface  waves  as  a  consequence  of  lateral 
variations.  Lateral  variations  in  structure  have  not,  of  course, 
been  taken  into  account  in  the  theoretical  calculations.  These 
effects  arc  quite  apparent  in  Figure  10,  however,  and  it  is  easy 
to  sec  that  in  the  vicinity  of  the  epicenter  the  pattern  most 
closely  approaches  the  theoretical  pattern,  while  at  greater 
distances  the  energy  has  been  refracted  and  scattered  by 
velocity  gradients  along  the  path  so  as  to  distort  the  basic 
symmetry  observed  near  the  source.  We  note  that  the  effects 
of  lateral  variations  in  structure  on  the  pattern  shapes  will 
be  frequency-dependent,  and  Figure  10  illustrates  this  fact. 

*Of  course  the  possibility  of  conversion  due  to  cracks  controlled 
by  the  local  stress  field  cannot  be  ruled  out. 
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One  of  the  major  considerations  of  this  study  is  the  com¬ 
parison  of  an  earthquake  field  with  those  from  the  explosions. 
The  Fallon  earthquake  was  located  very  near  to  the  SHOAL  epi¬ 
center  and  so  a  comparison  with  SHOAL  is  most  straightforward. 

Figure  12  shows  the  Love  wave  radiation  patterns  from  the 
Fallon  earthquake;  the  inset  shows  the  double  couple  equivalent 
which  gives  a  best  fit  to  these  observations.  The  same  equiva¬ 
lent  was  used  to  fit  the  Rayleigh  wave  patterns  and  the  L/R 
data  as  well,  of  course.  The  parameters  for  this  equivalent 
are  to  be  compared  to  the  double  couple  parameters  obtained 
from  the  SHOAL  explosion.  The  most  significant  difterence 
between  the  double  couple  equivalents  for  the  two  sources  is 
the  difference  in  depth  of  the  earthquake  double  couple.  In 
view  of  the  reasonably  strong  constraints  placed  on  depth  by 
the  radiation  patterns,  we  are  confident  that  an  appreciable 
difference  in  depth  exists.  We  therefore  have  an  cxanple  of 
what  a  triggered  event  at  depth  would  contribute  to  the  surface 
wave  radiation  from  an  explosion.  Here  we  note  that  the  earth¬ 
quake  has  a  significantly  lower  body  wave  magnitude  than  the 
explosion  (Table  I)  and  hence  is  energetically  of  the  size 
and  at  a  depth  that  has  been  hypothesized  for  triggered  events 
associated  with  explosions  in  the  range  of  SHOAL  and  BILBY. 

We  immediately  note,  however,  that  beside  the  fact  that  an 
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equivalent  double  couple  at  the  surface  fits  best  to  the  explo¬ 
sion  data  rather  than  one  at  significant  depth,  it  is  also 
evident  that  the  double  couple  for  the  earthquake  has  a  much 
higher  excitation  in  the  10  to  30  second  period  range  than 
does  that  for  the  explosion.  This  is  evident  by  comparion 
of  the  Love  wave  amplitudes  for  Fallon  with  those  for  SHOAL 
and  is  even  clearly  true  for  BILBY  which  is  more  than  an  order 
of  magnitude  larger  energetically  than  Fallon.  Thus,  we  would 
expect  that  if  a  triggered  earthquake  occurred  with  the  explo¬ 
sions  that  it  would  have  chracteristics  very  similar  to  the 
l-allon  earthquake,  especially,  of  course,  for  the  SHOAL  explo¬ 
sion.  However,  we  see  that  the  observed  characteristics  of  the 
anomalous  double  couple  components  associated  with  the  explo¬ 
sions  are  very  different  from  those  of  the  earthquake,  and 
conclude  in  accordance  with  Archambeau's  theory  that  the  most 
probable  mechanism  is  tectonic  release  associated  with  the 
highly  symmetric  fracture  zone  created  by  the  explosive  shock 
wave . 

One  other  important  conclusion  to  be  drawn  from  the 
difference  in  the  long-period  double  couple  excitation  of 
the  earthquake  relative  to  the  explosions  is  the  fact  that 
we  can  quite  easily  distinguish  between  earthquakes  and  explo¬ 
sions  just  on  the  basis  of  the  Love  wave  excitation  relative 
to  m^,  down  to  m^  -  4.4.  This,  of  course,  has  been  observed 
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before  (Lvernden,  1971),  at  least  indirectly,  in  terms  oi 

m.  vs.  M  differences.  While  it  is  clear  that  caution  is 
b  s 

necessary,  since  we  are  considering  only  one  earthquake  which 
may  have  somewhat  special  characteristics ,  it  seems  quite 
unlikely  that  these  differences  would  change  very  radically 
with  decreasing  magnitude.  Thus  we  expect  to  see  differences 
in  the  quadrupole  components  of  earthquakes  and  explosions, 
especially  in  the  Love  wave  excitation,  down  to  magnitudes 
below  four  for  most  if  not  all  earthquake-explosion  pairs.  We 
shall  discuss  this  point  again  in  the  next  section,  where  we 
will  consider  these  inherent  source  differences  in  terms  of 
Love  to  Rayleigh  wave  spectral  ratios. 

Another  aspect  of  the  comparison  between  the  quadrupole 
components  of  SHOAL  and  Fallon  is  the  similarity  in  the 
orientation  of  the  equivalent  double  couples.  We  expect  that 
if  the  origin  of  the  Love  waves  for  the  explosion  is  tectonic 
stress  relaxation,  then  the  orientation  of  the  double  couple 
equivalent  for  the  explosion  should  be  similar  to  that  for 
the  earthquake,  provided  of  course  that  the  prestress  field 
is  reasonably  well-behaved  and  such  that  averages  taken  over 
volumes  of  the  order  of  kilometers  in  diameter  do  not  vary 
much  in  this  region.  We  observe  from  comparison  of  tiie  para¬ 
meters  describing  the  double  couple  equivalents  for  the 
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explosion-earthquake  pair,  as  illustrated  by  the  insets  in 
Figures  10  and  12,  respectively,  that  there  is  a  strong 
correlation  of  the  double  couples  in  terms  of  their  orientations, 
although  there  is  a  large  difference  in  their  relative  magni¬ 
tudes,,  This  is  consistent  with  the  idea  that  tectonic  relaxa¬ 
tion  of  a  regional  stress  field  was  involved  in  both  cases, 
and  that  the  earthquake  was  much  more  efficient  at  low 
frequencies . 

I f  we  examine  the  radiation  patterns  shown  in  Figure  12 
in  more  detail,  we  see  that  the  theoretical  quadrupole  pattern 
shown  is  a  rather  rough  approximation  to  the  observations, 
although  a  good  first-order  fit.  There  are  at  least  two 
possible  reasons  for  this:  first,  we  expect  rupture 
propagation  effects  to  distort  the  patterns  at  wavelengths 
comparable  to  the  rupture  length  in  the  manner  described 
earlier;  second,  we  expect  structural  variations  to  cause 
perturbations  in  the  patterns.  Both  effects  appear  to  us 
to  be  present  here.  Note  in  particular  the  distortion  of 
the  pattern  at  T  =  9  seconds;  the  amplitudes  on  the  NE  lobe 
of  the  pattern  appear  to  be  3  to  4  times  as  large  as  those 
within  any  of  the  other  lobes.  This  effect  is  seen  to  diminish 
rapidly  with  increasing  period.  We  consider  this  deviation 
from  a  pure  quadrupole  pattern  to  be  predominantly  due  to 
rupture  propagation,  although  we  cannot  definitely  rule  out 
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structural  focusing  effects „  Although  the  amplitudes  were  low 
at  short  periods  for  the  SHOAL  event,  we  did  not  see  similar 
effects;  and  if  either  structure  or  fault  triggering  were 
involved  in  large  measure,  such  a  distortion  should  have 
been  evident„ 

Structural  effects  are  obviously  strong,  however,  and  we 
note  that,  as  with  the  observations  of  SHOAL  and  BILBY,  surface 
waves  propagate  efficiently  east  of  the  Rocky  Mountains,  It 
is  worth  noting  that  the  contour  pattern  swings  away  from 
north  toward  south  with  increasing  period  in  the  eastern  part 
of  the  continent  as  well.  This  implies  that  in  the  north  there 
is  a  higher  velocity  in  the  lower  crust  and  upper  mantle  compared 
to  the  south;  or  at  least,  such  a  variation  is  consistent  with 
this  sort  of  refraction  effect;  this  agrees  with  our  knowledge 
of  the  structural  variation  in  this  part  of  the  continent. 

Figure  13  shows  the  Rayleigh  wave  patterns  from  Fallon 
along  with  the  theoretical  pattern  produced  by  the  same  quadru- 
pole  used  to  fit  the  Love  waves.  In  this  case,  because  of  the 
orientation  of  the  equivalent  double  couple,  the  pattern 
shape  changes  with  frequency.  We  observe  a  rather  remarkable 
first-order  agreement  between  the  theoretical  and  observed 
patterns,  particularly  when  we  consider  the  combined  Love  and 
Rayleigh  wave  data.  Again  we  observe  strong  perturbations  in 
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the  patterns  due  to  structure,  especially  in  the  east.  However, 
close  to  the  source  we  again  observe  what  appears  to  be 
evidence  for  rupture  propagation  in  the  distortion  of  the 
pattern,  resulting  in  larger  amplitudes  to  the  north  at  short 
periods.  Again  a  strong  conclusion  to  this  effect  is  not 
possible,  but  at  least  the  Love  and  Rayleigh  waves  both  show 
the  effect,  as  would  be  expected  if  it  were  due  to  rupture 
propagation. 

In  the  east  the  effects  of  lateral  structural  variations 
are  different  for  the  Rayleigh  waves  than  for  Love  waves  in 
that  we  observe  no  similar  rotation  of  the  pattern  from  north 
to  south  for  Rayleigh  waves.  However,  the  Rayleigh  wave 
excitation  is  smaller,  and  noise  may  obscure  the  true  behavior 
of  the  pattern  as  a  function  of  period. 

A  composite  picture  of  the  fit  to  both  the  Love  and 
Rayleigh  wave  radiation  is  provided  by  the  fit  to  the  Love  to 
Rayleigh  spectral  ratio.  Figure  14  shows  the  L/R  radiation 
pattern  at  T  =  16  seconds,  a  period  of  high  amplitude  for 
both  types  of  surface  waves.  Considering  the  combined  effects 
of  lateral  structural  variations  and  rupture  propagation  as 
perturbing  influences,  the  simple  quadrupole  pattern  shown 
gives  a  remarkably  good  fit  to  the  data.  The  important  point 
here  is  that  the  L/R  ratio  for  Fallon  is  much  larger  than 
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the  L/R  ratio  for  SHOAL,  and  this  is  true  at  every  azimuth 
with  the  exception  of  only  those  (singular)  points  along 
the  node  lines,,  This  may  thus  potentially  be  a  useful  dis¬ 
criminant  between  explosions  and  earthquakes,  especially  If 
it  holds  over  a  wide  period  range.  In  the  next  section  we 
will  investigate  the  period  range  over  which  this  difference 
applies;  we  will  see  that  for  observations  at  larger  distances, 
it  is  necessary  to  consider  the  L/R  ratio  at  progressively 
lower  frequencies  in  order  clearly  to  distinguish  between  the 
two  source  types  on  the  basis  of  larger  L/R  for  earthquakes. 
Nevertheless  the  difference  holds. 

Spectral  comparisons 

In  the  previous  section  we  focused  most  of  our  attention 
on  the  azimuthal  variation  of  the  fields  from  these  sources, 
and  to  some  extent  on  their  distance  dependence.  In  this 
section  we  look  more  closely  at  the  spectral  properties  of 
the  files.  It  is,  of  course,  important  to  keep  in  mind  that 
the  spectral  characteristics  vary  with  azimuth  and  distance, 
and  we  will  therefore  attempt  to  integrate  the  results  of 
the  previous  section  with  the  present  discussion. 

Figure  15  shows  observed  surface  wave  spectra  from  BILBY 
at  progressively  larger  distances  along  a  line  of  stations  at 
roughly  the  same  azimuth.  The  radiation  patterns  for  BILBY  show 
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that  stations  along  this  azimuthal  sector  lie  near  the  maximums 
of  the  quadrupole  lobes  for  both  Love  and  Rayleigh  waves,,  These 
spectra  show  that  the  Rayleigh  wave  excitation  is  generally 
greater  than  the  Love  wave  excitation  over  the  range  of  obser¬ 
vations  from  5  to  40  seconds,  and  that  the  observed  spectra 
peak  in  a  range  near  10  seconds  period,,  The  exception  to  these 
generalizations  is  the  Love  wave  spectrum  at  BX-UT,  which 
shows  an  increase  in  amplitude  in  the  30  to  40  second  period 
range.  Since  the  time  sample  used  to  obtain  the  spectral 
estimate  at  this  near  distance  range  was  short  (the  signal 
being  quite  pulse-like)  ,  so  that  corrections  for  sample  trend 
and  the  like  are  therefore  more  uncertain,  we  feel  that 
this  feature  in  the  spectrum  may  be  an  artefact  of  our  procedure* 

The  BILBY  spectra  can  be  usefully  compared  with  the  Fallon 
earthquake  spectra  shown  in  Figure  16,  We  see  that  for  Fallon 
the  Love  wave  spectra  are  generally  greater  than  the  Rayleigh 
wave  spectra  at  all  distances,  which  is  opposite  to  the 
relation  for  the  explosion,  and  that  the  spectral  peak  is  in 
a  range  near  20  seconds,  although  the  Rayleigh  wave  spectrum 
does  not  have  a  very  clearly  defined  absolute  maximum.  In  any 
case,  even  a  casual  comparison  of  the  spectral  characteristics 
of  the  surface  wave  spectra  from  the  two  source  types  shows 
that  the  Love  wave  excitation  from  the  earthquake  is  relatively 
much  greater  than  from  the  explosion  and  that  the  spectral 
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excitation  is  shifted  toward  longer  periods.  Much  of  the 
detailed  character  of  the  spectra  shown  is  of  course  controlled 
by  the  structure;  however,  source  spectral  characteristics 
appear  to  result  in  changes  in  the  gross  properties  just 
mentioned . 

Figure  17  shows  the  comparable  spectra  for  the  SHOAL 
event.  In  general,  the  spectral  amplitudes  are  rather  low  and 
therefore  noise  is  a  significant  background  effect.  In  particu¬ 
lar,  we  regard  the  rise  in  the  spectra  beginning  around  30 
seconds  and  extending  to  the  40  second  limit  to  be  spurious 
and  to  be  due  either  to  imperfect  knowledge  of  the  instrument 
response  at  long  periods,  or  ambient  background  seismic  noise, 
or  both.  Ignoring  this  suspect  part  of  the  spectrum,  we  note 
that  the  spectra  of  both  types  of  surface  waves  peak  near  10 
seconds  (actually  at  periods  somewhat  less  than  10  seconds  as 
compared  to  BILBY,  where  the  spectral  maxima  occur,  on  the 
average,  at  periods  somewhat  greater  than  10  seconds).  The 
Rayleigh  wave  amplitudes  in  this  restricted  period  range  are 
of  the  order  of  or  greater  than  the  Love  wave  amplitudes.  Again 
we  observe  that  the  explosion  produces  surface  waves  with 
relatively  greater  excitation  at  short  periods  than  does  the 
comparable  earthquake,  but  does  not  produce  Love  waves  with 
anything  approaching  the  same  efficiency  as  a  comparable 
earthquake „ 
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A  more  quantitative  means  of  comparing  the  frequency 
dependence  of  the  Love  wave  excitation  between  an  earthquake 
and  an  explosion  is  to  calculate  the  ratio  of  the  normalized 
love  wave  spectrum  for  the  explosion  by  that  of  the  earthquake, 
figure  18  shows  the  Love  wave  spectra  for  SHOAL  and  Fallon 
recorded  over  nearly  identical  paths.  The  spectral  ratio  for 
the  Love  waves  should  therefore  be  essentially  independent 
of  path  and  should  reflect  source  differences  directly, 
furthermore,  the  Love  wave  excitation  is  a  very  slowly 
varying  function  of  depth  and,  for  a  source  depth  difference 
of  around  20  km  for  these  two  sources,  will  have  a  small 
effect  (Harkrider,  1964);  hence  we  see  principally  the  spec¬ 
tral  differences  in  the  source  equivalents  themselves.  The 
ratio  shown  in  Figure  18  indicates  the  relatively  greater 
short-period  excitation  of  Love  waves  for  an  explosion 
relative  to  an  earthquake  of  comparable  magnitude. 

In  terms  of  a  relaxational  source,  this  would  be  explained 
in  terms  of  a  much  larger  source  dimension  for  the  earthquake 
than  for  the  tectonic  source  associated  with  the  explosion. 

This  is  consistent  with  the  hypothesis  of  stress  relaxation 
from  the  region  around  the  shock- induced  shatter  zone.  We 
might  emphasize  that  even  though  the  P  wave  magnitude  of 
the  earthquake  was  significantly  smaller  than  the  explosion, 
the  source  dimension  appears  to  be  substantially  larger.  We 
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have  shown  a  1/T  curve  only  for  comparison  with  the  data; 
however,  we  observe  that  where  the  power  is  high  for  both 
terms  of  the  ratio,  the  data  appear  to  fit  such  a  dependence 
quite  well.  We  will  consider  theoretical  predictions  of  the 
relative  Love  wave  excitation  for  explosion  and  earthquake 
models  in  a  later  study,  with  the  objective  of  comparing 
these  predictions,  for  various  model  types,  to  data  such 
as  these. 

Since  the  power  at  the  long-period  end  of  the  Love  wave 
spectrum  for  the  SHOAL  explosion  was  low  and  probably  contami¬ 
nated  by  noise  to  a  significant  degree,  we  have  plotted  the 
sample  spectral  ratio  for  the  BILBY-Fallon  pair,  since  BILBY 
has  reasonably  high  power  in  the  20  to  40  second  period  range. 
Figure  19  shows  the  results,  with  a  reasonably  close  fit  to  a 
1/T  dependence  at  the  long-period  end  of  the  spectrum.  Taken 
together,  then,  the  spectral  ratios  for  Love  waves  from  these 
expl,  sions  relative  to  the  Fallon  earthquake  appear  to 
follow  a  dependence  which  is  roughly  1/T  in  the  range  which 
the  observations  are  reliable.  This  merely  quantifies,  somewhat, 
our  earlier  statment  that  explosion-generated  Love  waves  are 
predominantly  of  shorter  period  than  Love  waves  from  an  earth¬ 
quake  of  comparable  magnitude  (or  energy) .  An  explanation  is 
of  course  that  the  source  dimensions  are  quite  different,  the 
earthquake  being  the  larger. 
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Perhaps  the  most  interesting  result  of  this  study  is 
that  the  Love  wave  excitation  from  the  explosions  is  sig¬ 
nificantly  different  than  that  from  the  earthquake.  This 
difference  of  course  is  expected  for  an  idealized  explosive 
source;  the  only  question  is  how  strong  the  anomalous  effects 
might  be  in  actual  cases.  We  have  seen  that  while  the  anomalous 
field  appears  to  be  tectonic  in  origin,  it  is  not  so  over¬ 
whelming  as  to  mask  distinguishing  properties  of  the  purely 
explosion-generated  field. 

To  summarize  these  differences,  we  again  consider  the 
Love  to  Rayleigh  wave  ratio  (L/R)  for  these  events.  Figures  9 
and  11  show  the  spatial  character  of  this  ratio  for  the  two 
explosions  and  comparison  with  the  ratio  for  Fallon  in  Figure  14 
indicates  the  differences  in  these  characteristics  at  a  fixed 
frequency.  We  see  that  at  nearly  every  azimuth  the  L/R  ratio 
for  the  earthquake  is  large  compared  to  that  for  the  explo¬ 
sions  and  especially  so  along  the  strike  of  the  "fault  plane". 
This  suggests  that  by  checking  the  L/R  spectral  ratio  at  only 
a  few  stations  we  might  be  able  to  distinguish  between  earth¬ 
quakes  and  explosion  sources,  It  is  important,  however,  to 
investigate  the  detailed  behavior  of  this  ratio  as  a  function 
of  frequency  in  order  to  see  whether  the  difference  holds  over 
a  wide  range  of  frequency  and  whether  there  is  some  frequency 
where  the  difference  is  maximized.  In  particular,  does  L/R 
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for  earthquakes  remain  large  compared  to  that  for  explosions 
as  the  period  increases?  Also,  can  the  difference  in  the 
L/R  excitation  for  different  small  sources  be  distinguished 
at  large  distances? 

In  order  to  answer  these  questions  we  consider  the  L/R 
spectra  at  stations  along  a  nearly  constant  azimuth  to  the 
east  of  the  studied*  This  selection  of  azimuth  does  not 
maximize  the  difference  between  L/R  for  the  earthquake  compared 
to  the  explosions,  as  can  be  seen  from  again  comparing  Figures 
9,  11  and  14;  and  is  representative  of  the  average  difference 
likely  to  be  seen  over  a  distribution  of  stations  rather  than 
the  extreme  difference* 

Figure  20  shows  the  L/R  spectra  at  a  relatively  close 
range  for  the  three  events  studied*  It  is  clear  that  the 
closely  related  SHOAL-Fallon  pair  are  very  different,  since 
the  L/R  spectra  for  SHOAL  remains  rather  flat  and  near  unity 
(actually  less  than  unity),  while  that  for  Fallon  is  nearly 
three  at  the  10-second  limit,  increasing  with  increasing  period 
and  approximately  an  order  of  magnitude  larger  at  30  seconds 
period.  Generally  speaking,  L/R  for  BILBY  is  similar  to  that 
for  SHOAL,  although  on  the  average  it  is  somewhat  larger 
over  this  period  range.  Thus  over  the  frequency  range  where 
the  observations  are  reliable,  L/R  for  the  earthquake  is  much 
larger  than  for  the  explosions  and  increases  with  increasing 
period. 


-41- 


Figure  21  shows  the  L/R  spectra  for  the  events  at  an  inter 
mediate  teleseismic  range,  at  approximately  the  same  azimuths 
as  those  in  Figure  20.  We  note  that  roughly  the  same  relation¬ 
ship  between  the  events  again  holds.  In  this  case,  however, 
the  spectral  amplitudes  for  BILBY  appear  to  be  lower,  on  the 
average,  than  SHOAL.  However,  note  that  the  distinction  between 
the  earthquake  and  the  explosions  in  terms  of  L/R  at  the 
short-period  end  of  the  spectral  range  has  been  lost.  This  is 
probably  due  to  the  effects  of  structure,  e.g.,  scattering  from 
lateral  variations.  We  also  note  that  several  maxima  and  minima 
occur  in  the  spectra,  whereas  the  spectra  at  the  closer 
distance  were  more  regular.  Again  we  condsider  this  to  be 
associated  with  structural  irregularities.  This  emphasizes  the 
difficulties  involved  in  isolating  a  spectral  peak  (or  a 
minimum  in  the  spectrum  of  the  vertical  component  of  the 
Rayleigh  wave)  which  can  be  associated  with  a  source  depth 
effect.  Ilarkridcr  (1970)  considers  in  detail  the 
question  of  the  effect  of  source  depth  on  the  L/R  spectral 
ratio,  using  the  Fallon  earthquake  as  an  example. 

Finally,  Figure  22  shows  the  spectra  at  a  distance  of 
nearly  3000  km  from  the  sources.  We  note  that  the  spectrum 
from  Fallon  has  an  L/R  ratio  which  increases  dramatically 
with  period  and  actually  extends  to  a  value  considerably  larger 
than  20  at  the  longest  period  of  the  range.  However,  since  the 
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observed  power  at  the  extreme  end  of  the  range  was  low  for 
ballon,  we  do  not  consider  the  data  beyond  20  to  25  seconds 
to  be  very  reliable,,  Nevertheless  the  spectrum  maintains  a 
much  higher  value  than  the  explosions  over  a  considerable 
period  range,,  We  also  note  that  the  SHOAL  spectrum  exhibits 
a  large  peak  near  30  seconds  period;  however,  the  power  was 
extremely  low  at  this  distance  and  at  this  period,  and 
considering  the  noise  levels  and  the  effects  of  structure  on 
such  a  spectral  ratio,  we  conclude  that  this  peak  may  very 
well  be  spurious* 

Taken  together,  the  last  three  figures,  along  with  the 
radiation  patterns  shown  for  L/R  earlier,  show  that  the  L/R 
ratio  is  larger  for  this  earthquake  than  for  the  explosions 
over  a  period  range  from  around  10  to  35  seconds*  It  is  sig¬ 
nificant  that  this  holds  at  all  azimuths  (the  few  exceptions 
being  explained  by  nodes  in  the  radiation  patterns)  and  over 
the  distance  range  from  about  500  to  3000  km„  The  difference 
in  the  L/R  ratio  between  this  earthquake  and  the  comparable 
explosion  appears  to  increase  with  increasing  period  quite 
strongly,  ranging  from  a  factor  of  2  or  3  difference  near 
15  seconds  to  a  factor  of  around  10  near  30  seconds*  Various 
maxima  and  minima  appear  in  the  spectra,  depending  on  the 
distance  of  observation;  these  are  most  likely  associated 
with  scattering  due  to  structural  irregularities,  the  Love 
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waves  being  affected  differently  than  the  Rayleigh  waves. 

This  strongly  suggests  that  determination  of  depth  from 
spectral  minima  (or  maxima)  is  probably  decidedly  difficult 
and  perhaps  not  useful  in  a  practical  sense.  This  conclu¬ 
sion  is  consistent  with  that  obtained  by  Harkrider  (1970)  in 
his  detailed  study  of  the  surface  wave  spectra  and  L/R 
ratios  from  Fallon. 

We  conclude  that  the  L/R  spectral  ratio  may  be  a  reliable 
discriminant  for  magnitudes  near  4.5,  and  possibly  to 
lower  magnitudes.  We  have  not,  of  course,  demonstrated 
that  this  difference  is  appropriate  for  all  earthquakes,  and 
it  is  entirely  possible  that  some  earthquakes  may  have 
properties  quite  different  from  the  one  studied,  and  thus 
produce  L/R  spectra  similar  to  that  from  explosions.  In  this 
regard,  deep  earthquakes,  earthquakes  of  volcanic  origin,  and 
dip-slip  faulting  will  probably  not  show  such  a  strong  L/R 
difference  from  explosions  selected  from  arbitrary  tectonic 
regions.  However,  it  seems  clear  that  for  an  explosion  in  a 
tectonic  environment  which  produces  dip-slip  faulting,  we 
would  expect  a  corresponding  and  proportionate  decrease  in  Love 
wave  excitation  (since  such  excitation  is  also  due  to  a 
tectonic  origin)  as  the  earthquake,  and  hence  for  other  events, 
within  the  same  tectonic  region  an  L/R  difference  between 
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explosions  and  earthquakes  similar  to  that  observed  here  for  a 
region  typified  by  strike-slip  faulting.  In  any  case,  the  use 
of  the  L/R  spectral  ratio  as  a  discriminant  appears  potentially 
useful,  but  the  possibility  of  practical  application  needs 
additional  study,  using  a  wider  selection  of  earthquakes  and 
explosions  from  different  tectonic  environments  than  that  used 
here. 

A  second  conclusion  which  can  be  drawn  from  the  nature  of 
the  L/R  spectra  verifies  that  inferred  from  the  comparison  of 
only  the  Love  wave  spectra  for  the  events  as  given  in 
Figures  18  and  19.  In  particular,  the  in  ease  of  L/R  with 
respect  to  period  for  the  earthquake,  and  the  generally  higher 
value  of  this  ratio  over  the  entire  period  range  compared  to 
the  explosions,  indicates  that  the  anomalous  shear  wave  radia¬ 
tion  from  the  explosion  is  consistent  with  a  source  mechanism 
which  is  of  much  smaller  dimension  than  that  of  the  earthquake. 
It  is  important  to  note  here  that  the  effect  of  the  difference 
in  source  depth  between  the  earthquake  and  the  explosion  (20  km) 
is  such  as  to  predict  that  the  L/R  ratio  will  have  its  largest 
difference  at  short  periods,  with  the  shallow  source  excita¬ 
tion  being  the  greater,  and  that  the  L/R  ratios  become  nearly 
the  same  at  longer  periods  (30  seconds).  (See  Toksoz  et  al., 
1965,  Figure  6  and  13).  This  is  just  the  opposite  of  what  is 
observed,  and  hence  the  observation  is  probably  due  to  an 
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intrinsic  source  property  rather  than  a  medium  response  effect 
due  to  depth  of  source,  if  we  assume  that  the  origin  of  the 
anomalous  field,  which  manifests  itself  principally  as  l.ovc 
wave  radiation,  is  of  tectonic  origin,  then  we  arc  led  to 
conclude  that  the  effect  is  due  to  stress  relaxation  around 
the  shock  created  shatter  zone,  in  view  of  the  implied  small 
dimensions  of  this  source. 
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CONCLUSIONS 

Wo  draw  the  following  conclusions  from  this  study: 

(1)  The  earthquake  radiation  field  is  to  first  order  a 
quadrupolo  field  in  the  spectral  range  where  typical  wave¬ 
lengths  of  the  radiated  field  are  larger  than  the  rupture 
dimensions.  At  shorter  periods  there  are  indications,  in  hath 
the  radiation  patterns  and  spectra,  of  higher  order  multipole 
contributions  produced  by  rupture  propagation.  These  results 
arc  in  quantitative  agreement  with  the  theoretical  predictions 

(2)  Tectonic  stress  relaxation  is  by  far  the  most  likely 
origin  of  the  anomalous  contribution  to  the  radiation  field 
from  the  explosions.  From  our  analysis  of  the  spectral  and 
radiation  pattern  observations,  we  conclude  that  the  effect 
is  most  likely  due  to  stress  relaxation  aroung  the  shock- 
generated  shatter  zone  for  the  events  studied. 

(3)  The  anomalor.s  field  from  the  explosion  is  purely  quadrn- 
polc  over  a  wide  spectral  range.  Since  the  Rayleigh  wave  radia¬ 
tion  patterns  do  not  change  toward  being  a  qudrupolc  type  at 
longer  periods,  and  since  the  Love  to  Rayleigh  excitation  ratio 
is  near  unity,  we  conclude  that  the  source  spectrum  for  the 
anomalous  field  is  similar  in  shape  to  that  for  the  explosive 
monopole  component  of  the  total  field.  This  is  in  agreement 

wit)  theoretical  predictions  for  stress  relaxation  around  the 
shock-generated  failure  zone. 
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(4)  The  spectral  ratios  of  Love  waves  from  the  explosions, 
relative  to  the  earthquake,  show  u  long-period  excitation  for 
the  anomalous  explosion  field  which  has  a  fundamentally  different 
frequency  dependence  within  the  range  of  observation  (10  to  40 
seconds  period):  the  anomalous  explosion  field  has  a  frequency 
dependence  resembling  um+*  (where  m  is  some  positive  integer), 
while  the  earthquake  spectrum  has  a  dependence  like  um.  Thus 

t he  ratio  decreases  us  u  with  increasing  period  in  this  range. 
This  implies  that  the  source  of  the  anomalous  explosion  field 
is  quite  different  from  the  comparable  earthquake,  and  we 
conclude  that  if  the  origin  is  tectonic,  it  is  associated  with 
a  region  of  relatively  small  dimension.  This  suggests  that 
these  effects  are  due  to  stress  relaxation  around  the  shatter 
tone. 

(5)  The  ratio  of  Love  to  Rayleigh  wave  spectral  amplitudes 
generated  by  the  sources  studied  here  show  that  the  eurthquake 
(m^  =  4.4)  excited  Love  waves  much  more  efficiently  than  did 
the  explosions,  with  the  difference  in  excitation  increasing 
with  increasing  period.  We  conclude  that  this  difference  could 
be  useful  in  discrimination  of  the  two  source  types,  even  for 
relatively  low  magnitudes,  e0g„,  near  m^  =  4.4,  We  also  con¬ 
clude  that  the  process  of  tectonic  release  was  different  for 
the  explosions  than  for  the  earthquake,  which  was  of  comparable 
body  wave  magnitude.  This  strengthens  our  prior  conclusion 
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that  any  tectonic  effects  are  due  to  stress  relaxation 
around  the  explosion-generated  failure  zone. 
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Figure  1.  Distribution  of  stations  used  in  the  analysis  of  the  SHOAL,  Fallon  and  BILBY  events 


Figure  2.  Long- 


Group  velocity, 
km/sec 


Figure  3.  Staples  of  surface  wave  signals  and  the  associated 
group  velocity  dispersion  obtained  froa  these 
signals  using  narrow  band  filtering  Methods. 


figure  4.  samples  of  surface  wave  spectra  fro*  the  three  events  studied.  The  signal  sample  and 
the  group  velocity  for  the  BUSY  Rayleigh  spectra  arc  shown  in  Figure  3.  Fourier  spectra  of  the 
Ignats  are  shown  on  the  left  and  the  group  spectra  obtained  by  narrow  filtering  on  the  right. 


8  *  Dip  angle  0  *  Azimuth  angle 

4>  *  Strike  angle  n  *  Surface  normal 

X  «  Slip  angle  m,/«  Orthogonal  vectors 

m  plane  of  surfoce 

t^°*  «  Traction  vector  in 
plane  of  surfoce 

Figure  S.  Fault  (or  traction)  plana  geometry.  The  point  (0)  * 

it  an  arbitrary  point  of  observation,  while  the 
vector  K  is  in  the  direction  of  geographic  North. 


Figure  6,  Variation  of  the  Love  to  Rayleigh  (L/R)  wave  ratio 
aa  a  function  of  X,  the  (double  couple)  "slip 
angle",  for  superposed  aonopole  and  quadrupolc 
(double  couple)  coaponents.  The  relative  excitation 
of  the  quadrupolc  to  aonopole  coaponent  is  fixed 
at  F  •  O.SI  for  all  frequencies. 


THEORETICAL  LOVE  TO  RAYLEIGH 
WAVE  AMPLITUDE  RATIOS  FOR  AN 
EXPLOSION  AND  DOUBLE  COUPLE 
AT  0  km  DEPTH. 


Figure  7.  Variation  of  the  Love  to  Rayleigh  (L/R)  wave  ratio 
as  a  function  of  {,  the  (double  couple)  "dip  angle", 
for  superposed  nonopole  and  (double  couple)  ouadrupole 
coaponents.  All  other  source  paraneters  are  fixed,  as 
indicated. 
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Figure  9.  The  theoretical  Love  to  Rayleigh  (L/R)  wave 
ratio  for  BILBY,  at  T  •  IS  seconds,  as  a 
fuention  of  aziauth  using  the  source  paraaeters 
of  Figure  B.  The  observed  ratios  at  a  nuaber  of 
stations  with  good  signal  to  noise  ratios  for 
both  Love  and  Rayleigh  waves  are  shown  at  their 
appropriate  aziauths.  This  also  shows  the  nature 
of  the  L/R  aziauth  variation  and  that  a  single 
equivalent  source  fits  both  Love  and  Rayleigh 
wave  radiation  for  this  event. 
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Figure  il.  The  theoretical  Love  to  Rayleigh  (L/R)  wave 
ratio  for  SHOAL,  at  T  ■  15  seconds,  as  a 
function  of  azimuth  using  the  source  parameters 
of  Figure  10.  Observations  at  various  azimuths 
where  signal  to  noise  ratios  were  high  are 
indicated  and  identified  by  abbreviated  station 
symbols.  This  also  shows  the  nature  of  the  L/R 
azimuth  variation  (which  is  nearly  identical  to 
that  for  BILBY)  and  that  asingle  equivalent 
source  fits  both  Love  and  Rayleigh  wave  radiation 
for  this  event. 
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Figure  14.  The  theoretical  Love  to  Rayleigh  (L/R)  wave 
ratio  for  Fallon,  at  T  *  16  seconds,  as  a 
function  of  azimuth  using  the  source  parameters 
obtained  for  the  results  shown  in  Figures  12 
and  13.  Observations  at  various  azimuths  where 
the  signal  to  noise  ratios  were  high  are  indi¬ 
cated  and  identified  by  abbreviated  station 
symbols.  This  shows  the  nature  of  the  L/R 
azimuth  variation  for  an  earthquake  in  the  sane 
environment  a:;  the  comparable  explosion:  SHOAL 
(see  Figure  11).  The  L/R  ratio  for  this  t'vent  i  *  much 
larger  than  that  for  the  explosion  at  all  com¬ 
parable  azimuths,  except  along  nodal  lines. 


Amplitude,  /x/cps 


Period,  sec 


Hire  1  .  Surface  wave  spectra  at  3  distances  for  the 
KILBY  event.  The  points  of  observetlon  are  at 
approximately  the  sane  azimuth  from  the  soyree. 
Note  that  the  Rayleigh  wave  excitation  is 
generally  higher  than  the  Love  wave  excitation. 
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Figure  li.  Surface  war#  spectra  at  S  distances  for  the  Fallon 

event.  The  Mints  of  observation  are  at  approxinately 
the  sane  azlnuth  free  the  source.  Note  that  the  Love 
wave  excitation  is  general 1/  greater  than  the 
Kajrleigh  wave  excitation. 
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figure  1 7*  Surface  wave  spectra  at  3  distances  for  the  SHOAL 
event.  The  points  of  observation  are  at  approxi- 
nately  the  sane  asinuth  fron  the  source.  The  Love 
and  Rayleigh  excitation  Is  of  roughly  the  sane 
order,  although  at  the  nearest  (and  no si  reliable) 
station  the  Rayleigh  wave  excitation  Is  generally 
greater  over  the  reliable  period  range. 
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Figure  It.  Love  wave  spectra  and  spectral  ratios  for  the 
SHOAL  and  Fallon  events.  The  upper  two  graphs 
show  the  Love  wave  spectra  for  the  individual 
•vents  with  the  Rayleigh  wave  spectra  given  for 
comparison.  Roth  spectra  are  individually 
normalised  to  unity  at  the  maximum  power  point 
ant  the  resulting  normalised  ratio  of  MIOAL  to 
Fallon  Love  wave  spectral  density  i*  plotted. 

The  sanpled  spectral  ratio  is  shown  by  the  data 
points,  the  dotted  curve  is  a  smoothed  fit  to 
these  sample  points  and  the  solid  curve  a  I/T  curve 
given  for  comparison. 
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I ifurc  19.  Love  wave  spectra  and  spactral  ratios  for  the 
•ILftY  aid  Fallon  events.  Tk<  upper  too  graphs 
shoo  the  Loro  wave  spectra  for  the  two  events 
and  the  lower  graph  the  normalised  spectral 
ratio  as  described  in  Figure  11. 


Figure  20. The  variation  of  L/R  versus  period  for  the  RILftY 
and  SHOAL  explosions  and  the  Fallon  earthquake, 
in  the  distance  range  near  <00  kn.  The  L/R  for 
the  earthquake  is  larger  than  2  and  increases 
with  increasing  period,  while  L/R  for  the  earth* 
quake  is  larger  than  2  and  increases  with  increas¬ 
ing  period,  while  L/R  for  the  explosions  is 
between  1  and  2  and  reaains  nearly  constant  with 
increasing  period.  This  difference  in  behavior 
is  a  consequence  of  the  nuch  larger  source 
dinension  associated  with  the  earthquake. 


rigure  21,  The  variation  of  l/t  versus  period  for  tha  BILBY 
and  SHOAL  explosions  and  the  Fallon  earthquake, 
in  tha  diatanca  range  near  1600  ka,  Tha  L/R  for 
tha  earthquake  it  again  considerably  larger  than 
that  for  tha  explosions  and  increases  with 
increasing  period  at  this  larger  distance  range, 
as  expected.  Kota,  however,  that  the  details  in 
the  variation  of  L/R  period,  i.e.,  local  naxina 
and  ninina,  are  considerably  different  than  in 
Figure  20  and  that  the  sepa  ation  in  the  L/R  is 
obscured  at  the  short  period  end  of  the  range. 
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Figure  22.  The  variation  of  L/R  versus  period  for  the  BILBY 
and  SHOAL  explosions  and  the  Fallon  earthquake , 
in  the  distance  range  near  2800  km.  At  large 
distances  the  L/R  separation  between  the  earth¬ 
quake  and  the  explosions  is  maintained  at  long 
periods.  The  peak  in  the  spectral  ratio  for 
Shoal  at  30  seconds  is  probably  due  to  a  path 
effect. 


APPENDIX  I 


T1  IliORliT I  CAL  BACKGROUND 


We  begin  with  the  stress  relaxation  source  described  by 
Archambeau  (1968)  and  Archambeau  and  Minster  (1971)  „  l  irst 
we  present  briefly  the  analytical  nature  of  this  field  repre¬ 
sentation,  and  show  that  for  the  purposes  of  the  present  study 
it  is  sufficient  to  replace  the  actual  source  by  an  equivalent 
point  source  whose  primary  contribution  to  the  far-field  radia¬ 
tion  is  a  quadrupole  term  at  the  long-wavelength  end  of  the 
spectrum.  We  are  thus  able  in  this  approximation  to  make  use  of 
the  point-source  models  described  by  Ben-Menahem  and  llarkrider 
(1964),  and  in  particular  to  use  llarkrider's  computer  program 
to  calculate  a  theoretical  fit  to  our  observations  of  radia¬ 
tion  patterns. 

For  the  explosion  sources  studied  here,  we  adopt  the 
hypothesis  that  the  observed  anomalous  shear  wave  energy  is 
due  to  tectonic  effects.  This  idea  has  been  considered  in  some 
detail  by  Press  and  Archambeau  (1962),  Toksoz  et  al„  (1965), 
Archambeau  and  Sammis  (1970),  and  Tsai  and  Aki  (1971),  among 
others.  At  the  present  stage  we  do  not  assume  any  specific 
model  for  the  process  of  tectonic  release,  i.e.,  we  do  not 
choose  beforehand  among  the  following  possibilities:  (1)  all 
the  tectonic  release  is  due  to  the  creation  of  a  simple, 
nearly  spherical,  shatter  or  failure  zone  caused  by  the  shock 
wave  from  the  explosion;  (2)  all  or  most  of  the  tectonic  energy 


release  is  due  to  the  triggering  of  rupturing  which  results  in 
the  creation  of  lineal  failure  zones  larger  than  the  shatter 
zone.  Instead  of  adopting  either  of  these  hypotheses,  we 
present  here  a  theory  which  covers  both  cases,  and  study  the 
observed  spectral  data  from  the  explosions  in  order  to 
attempt  to  determine  whether  or  not  tectonic  release  actually 
contributes  to  the  radiation  field,  and  if  so,  which  of  these 
"geometrical"  mechanisms  _s  responsible. 

As  far  as  the  direct  contribution  to  the  radiation 
field  fi"om  explosion-generated  shock  waves  is  concerned,  it 
will  be  assumed  as  a  working  hypothesis  that  the  field  is 
adequately  represented  by  a  monopolc  contribution.  Because 
of  inhoinogencit  ics  and  anisotropy  in  the  source  region, 
which  will  affect  the  conversion  of  the  shock  waves  to  elastic 
waves,  there  will  undoubtedly  he  higher-order  multipolc  contri¬ 
butions,  but  it  is  reasonable  to  consider  these  as  second  order 
compared  to  the  monopole  term  (Toksoz  et  al.,  1965),  and  also 
of  second  order  compared  to  the  hypothesized  tectonic  effects. 
However,  one  of  the  purposes  of  this  study  is  to  assess  the 
relative  importance  of  any  tectonic  effects  compared  to  the 
non-monopolc  contributions  to  the  far-ficld  radiation  from 
the  explosive  sources.  Consequently,  although  we  adopt  a 
reasonable  working  hypothesis  in  order  to  have  a  basis  for 
the  development  of  a  theoretical  framework,  we  do  not  exclude 


other  possibilities  from  our  interpretive  considcrut ions. 

Stress  relaxation  effects  for  both  explosions  and  earth* 
quakes  have  been  considered  in  some  detail  by  Archamheau 
(1968,  1971)  and  here  we  need  merely  summarize  these  briefly. 
The  treatment  is  based  on  an  ini t ial-value  formulation,  w  th 
the  medium  assumed  to  be  in  a  prestressed  equilibrium  stale, 
so  that  upon  creation  of  a  failure  lone  the  solid  medium 
relaxes  to  a  new  equilibrium  state  and  a  new  state  of  stress. 
This  is  accomplished  by  radiation  of  energy  away  from  the 
volume  surrounding  the  failure  zone,  with  a  net  reduction  of 
potential  (or  strain)  energy  within  this  region;  this  process 
is  naturally  described  as  an  initial-value  problem.  The  energy 
release  is,  of  course,  due  to  the  change  in  effective  material 
properties  in  the  failure  zone,  since  afterward  the  medium 
will  not  be  able  to  sustain  the  same  stress,  i.c.,  it  cannot 
maintain  the  initial  equilibrium. 

We  can  treat  such  a  problem  quite  simpiy  in  an  unbounded 
medium,  and  then  add  to  the  result  the  solution  to  the  homo¬ 
geneous  problem  (no  source  terms  or  initial  values)  in  ordev 
to  satisfy  specified  conditions  on  boundaries  within  or 
surrounding  the  actual  volume  of  interest. 

The  equations  of  motion  for  the  infinite  homogeneous, 
isotropic,  perfectly  clastic,  prestressed  solid  medium, 
neglecting  hiohor-ordor  interactions  with  the  prestress  field. 


y-sy 


can  he  reduced  to  the  set  of  wave  equation*  (Archanbeau,  I9?l): 


l  •  •  i.  2.  2.  4 


(1) 


where  the  four  potential*  representing  the  three  component* 
of  rotation  and  the  »ca!ar  dilatation  (the  scalar  having  index 
a  9  41)  provide  the  displacement  field  via  the  relation*: 


i}  , 
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ijk  \ 


(2) 


where  the  indice*  i,  j,  l  run  over  the  range  1,  2,  3  and 

i*  the  alternating  tensor.  Here  we  employ  the  usual  summa¬ 
tion  convention:  i»  the  shear  velocity  for  each  k  *  1,  2,  3, 

while  Vj  is  the  compressional  velocity,  q  ^  in  equation  (1) 
represents  any  explicit  body  force  (e.g.,  gravity  or  internally 
applied  force  loads).  In  the  present  problem  we  can  represent 
the  explosion  itself  by  a  distribution  of  forces  over  a  surface 
within  the  medium,  and  we  could  use  a  delta  function  form  for 
q^  to  acconplish  this  representation. 

The  potentials  arc  dynamic  field  potentials  defined 
(or  measured)  relative  to  the  final  equilibrium  state  of  the 
medium.  The  prestress  does  not  appear  explicitly  In  these 
equations,  but  it  will  appear  as  a  generalised  Initial  value 
for  these  potentials  since  they  have  non-sero  values  appropriate 


fl-5 


to  the  Initial  state  of  stress  at  the  onset  of  the  failure 
process  and  util  eventually  assume  final  (rcro)  values  appro¬ 
priate  to  a  state  of  stress  In  equilibrium  with  the  completed 
failure  tone.  The  chance  in  the  potentials  at  a  given  stage 
of  the  rupture  tone  configuration  is  given  by  Archanhcau  ( I9“l ) 
as: 


vv  •  o 

*•» 

and  thus  the  change  in  the  >  for  a  change  in  naterinl  proper- 
tics  due  to  failure  Is  always  a  harnonic  function. 

A  Green's  function  solution  for  equation  (I)  provides  a 
convenient  Method  of  introducing  the  initial  values  of  >  t 
and  shows  the  effects  of  other  source  terns,  as  well  as  inter¬ 
actions  between  the  field  and  boundaries  within  the  nedlun 
he  have  (c.g.,  Archanhcau,  1971): 

<1.  O  •  4tm  J  Co<£,  t;  r^,  t0)  ^  (r^.  t#>  4 ^ 

o  * 


vhcrc,  as  previously  noted,  we  will  use  the  Green's  function 
for  an  Infinite  space,  given  bv: 

C.  -  7T  Ifr  6.  •  4t  4<X-I.>  4(t*to> 

ft  9 

The  first  tern  in  (4)  corresponds  to  the  particular 
effects  of  body  force  terns  or  tine-varying  applied  forces, 
and  in  considerations  of  tectonic  effects  per  se  these  will 
he  onlued  iron  further  attention.  The  second  tern  will  also 
be  ignored  for  the  present,  since  it  represents  the  Inter¬ 
action  of  the  field  with  boundaries,  including  the  boundary 
enclosing  the  growing  failure  tone.  Since  the  individual  parts 
of  the  field,  both  due  to  the  source  and  to  boundary  "scattering" 
effects,  are  additive,  we  can  superpose  the  effects  of  the 
boundaries  later,  using  a  solution  involving  an  eigenfunction 
expansion  which  Is  equivalent  to  the  surface  Integral  tern, 
he  are  left  only  with  the  last  tern  In  the  solution  (4)  which 
is  thus  seen  to  take  into  account  the  Initial-value  effects. 

Now  we  consider  the  effect  of  introducing  a  shock™ 
generated  failure  rone  due  to  an  explosion  (for  details  of  the 
developncnt,  see  Archanheau,  1971).  he  assume  that  the  creation 
of  this  rone  proceeds  at  a  rate  faster  that  any  of  the  intrinsic 
wave  velocities  within  the  nediun,  and  if  we  take  the  shape  of 
the  rone  of  failure  to  be  spherical  as  a  first  approx inat ion, 


then  we  have  in  effect,  an  instaneou*  creation  of  a  shatter 
zone  In  thii  ca>e,  denoting  the  last  tern  in  (4)  »>  ,fl  r,  t|, 
we  have: 
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with  an  initial  Mate  given  l»> : 
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where  ft  is  the  ridiu*  of  the  shatter  rone  and  vJ?  is  the 
rupture  velocity,  which  we  here  take  to  he  equal  to  the 
shock  velocity  for  the  nediun.  Thus  the  initial  values  art 
defined  at  the  tine  of  conpiction  of  the  rupture,  w.id  re  I  is  i- 
proceeds  fron  that  line  onward  at  a  rate  governed  by  the 
intrinsic  velocities  within  the  nediun.  he  note  that  the 
initial  value  x*  is  Riven  by  the  solution  of  (J)  along  with 
continuity  conditions  on  the  rupture  surface  and  a  condition 
requiring  i*  to  approach  the  unperturbed  initial  value  at 
large  distances  fron  the  rupture  rone.  That  is,  wc  require 
a  solution  to  the  standard  prohlen  for  a  spherical  inclusion 
in  a  stressed  nediun.  Integrating  over  the  source  tire  variable 


t  *nJ  taking  account  of  the  discont inuous  nature  of  the  integrand. 


vo  find: 


(S) 


Here  *•  i s  a  known  harmonic  function  and  (.  is  the  Green's 

*  f$ 

function  for  the  inf  ini te-spac#  prohlen  Thus  we  have  an 
explicit  volume  intcitrai  to  evaluate,  which  is  a  quite 
straightforward  prohlen. 

Hhen  the  rupture  velocity  is  less  than  the  intrinsic 
velocity  v  ,  relaxation  of  stress  occurs  as  the  rupture 
volume  increases  or  changes  shape.  This  is  the  case  for 
spontaneous  failure  and  thus  is  appropriate  for  the  earth¬ 
quake  source  nodel. 

If  we  view  this  process  as  a  sequence  of  instananeous 
suh-proccsscs,  each  separated  by  .an  increnental  tine  5 tQ , 
we  can  add  together  solutions  of  the  forn  ($)  to  get  the 
effect  of  finite  rupture  velocity  in  the  stressed  neditin. 
In  the  Unit  as  't  approaches  zero,  this  will  ylc»ld  the 
required  solution  (sec  Archanbeau ,  1961  and  1971,  for  .an 
expanded  discussion),  Thus,  taking  the  rupture  growth  in 


isich  (towth  invrmnt  l<tin, 


step*  at  t  ^  lines  l  •  1,  I,  .  , 
sep#*1  »tcd  in  line  In-  we  obtain  by  tupnpo«ine  xoluti  tt, 
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llcrc  we  have  inserted  4t  cxpli.lt  ly  by  multiplying  at  I 

dividing  (S)  by  the  tine  increment  The  quint  it  *  lx  %nn 
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and  this  corresponds  to  the  change  in  the  cqul  I  ibr  lur,  field 
due  to  the  Increnrni  of  rupture  growth  which  occurs  at  t Ini 

lok* 

Sow  if  we  allow  &t(  to  approach  zero  while  *  I  mil t ancon* I y 

allowing  the  number  of  growth  increments  to  approach  infinite, 

the  summation  becomes  an  integral  and  the  difference  rati 

4»*/4t„  becomes  a  partial  derivative,  and  so  we  have: 
no  1 
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Here  H(t0»t0)  is  a  unit  step  function  and  :  is  the  tine 
during  which  rupturing  occurs,  Mien  tQ  >  t  this  step  function 


vanishes;  it  l*  unity  for  tQ  <  anJ  thus  it  nercly  serve* 
to  define  the  fine  interval  of  the  rupture  process  and  hence 
the  iutt rial  of  integration  over  the  source  ti»c  tQ. 

In  both  the  forns  (S)  and  (b),  the  solution  for  the  field 
cause/  h)  the  relaxation,  •  *,  has  the  sane  neaning.  In  particular, 
it  is  the  equilihriun  field  for  the  prestressed  nediun  with  an 
inclusion  (the  rupture  zone).  The  nature  of  the  inclusion,  l.e, 
its  shape,  size,  and  so  forth,  is  dependent  on  the  "source" 
tine  tQ,  since  we  allow  the  inclusion  to  change  shape  as  a 
function  of  tine  t  ,  Thus  >*  is  the  equilihriun  field  the 
nediun  would  relax  toward  if  the  rupture  naintained  the 
specific  forn  it  had  at  the  tine  t0> 

The  basic  integral  representations  in  (S)  and  (6)  can  he 
evaluated  in  a  relatively  straightforward  fashion  for  the  case 
in  wKch  the  naierial  is  unifomly  prestressed  in  shear  and 
when  the  basic  failure  has  a  spherical  shape.  Thus,  for  either 
instantaneous  spherical  shatter  zone  growth  or  for  rupture 
growth  at  a  rate  less  than  the  shear  velocity,  we  can  obtain 
closed-fora  solutions  (sec  Archanbeau,  1961  and  1971).  further, 
we  wan  allow  the  spherical  zone  to  grow  and  to  translate  In 
tine  within  the  stressed  nediun,  thereby  sinulating  a  propagating 
failure  zone  which  heals  itself,  l.e.,  locks  up  a  short  tine 
after  the  initiation  of  the  failure.  Thus  we  can  calculate  the 
field  fron  a  transient  type  of  failure  which  sweeps  through  a 


large  region  (sec  lor  details  Archanhcau  and  Minster,  1D71J  All 
these  solutions  have  the  same  analytic  form;  regarded  as  multi¬ 
pole  fields  in  the  frequency  domain  they  arc  given  by: 
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(see  Archambeau,  1968  and  1971).  The  quantities  are  the 

constant  components  of  the  pure  shear  prestress  field  we 
assume,  and  the  matrices  are  indexed  on  a  =  1,  2,  3,  4,  and 
m  =  0,  1,  2.  The  prestress  matrices  and  the  functions 
and  11^  determine  the  nature  of  the  multipole  coefficients 
and  .  These  functions  are  dependent  on  the  type  of 
rupture  involved,  and  in  general  are  rather  complicated 
series  of  hypergeometric  functions  (derived  in  detail  by 
Archambeau,  1968  and  1971,  and  Archambeau  and  Minster,  1971). 


